arXiv:1508.03588vl [cond-mat.soft] 14 Aug 2015 


Surface Tension regularizes the Crack Singularity of Adhesive Contacts 


Stefan Karpitschka^, Leen van Wijngaarden^, and Jacco H. Snoeijer^’^ 

^Physics of Fluids Group, Faculty of Science and Technology, 

Mesa+ Institute, University of Twente, 7500 AE Enschede, The Netherlands. 

^Mesoscopic Transport Phenomena, Eindhoven University of Technology, 

Den Dolech 2, 5612 AZ Eindhoven, The Netherlands 
(Dated: August 17, 2015) 

The elastic and adhesive properties of a solid surface can be quantified by indenting it with a 
rigid sphere. Indentation tests are classically described by the JKR-law when the solid is very stiff, 
while recent work highlights the importance of surface tension for exceedingly soft materials. Here 
we show that surface tension plays a crucial role even for stiff solids: it regularizes the crack-like 
singularity at the edge of the contact. We find that the edge region exhibits a universal, self-similar 
structure that emerges from the balance of surface tension and elasticity. The similarity theory 
provides a complete description for adhesive contacts, reconciling the global adhesion laws and local 
contact mechanics. 


The adhesion between two solid bodies in contact is 
extremely common in nature and technology [IHl]. Ad¬ 
hesive contacts are described by a classical law derived 
by Johnson, Kendall and Roberts (JKR) [3, providing 
the benchmark to characterize elastic and adhesive ma¬ 
terial properties EHH]. Despite its success, JKR theory 
does not provide a complete description of the physics 
inside the contact. Namely, the elastic problem is con¬ 
sidered without explicitly treating adhesive interactions, 
while furthermore the contact exhibits a crack-like sin¬ 
gularity at edge. This issue was elegantly resolved in a 
macroscopic theory, where, in analogy to fracture me¬ 
chanics, the elastic energy released by opening the crack 
is balanced by the work of adhesion [MS]. 

Recent studies on very soft gels and rubbers provided a 
very new perspective on adhesive contacts [HUTU] . Just 
like fluids, these soft materials are highly susceptible to 
a Laplace pressure due to surface tension [20H2l]- This 
insight exposed a profound link between “solid adhesion” 
and “liquid wetting” [IlHIS]: both are adhesive contacts, 
but whether they are solid-like or liquid-like depends on 
the contact size {£ in Fig.[^ compared to the elastocap- 
illary length 7 //i, the ratio of surface tension 7 to shear 
modulus /i. Liquid-like contacts were for example found 
for soft gels m and nanoparticles m- 

Intriguingly, the limit of liquid wetting does not suffer 
from the crack singularity: it is governed by a benign 
boundary condition, the wetting angle 9. This stark con¬ 
trast to the JKR-singularity, in combination with recent 
experiments |25j , once more provokes the question: What 
happens at the edge of adhesive contacts? 

In this Letter we show that surface tension is crucial 
also for stiff, adhesive contacts. Our analysis reveals how 
surface tension regularizes the crack singularity on a scale 
7 //i, replacing it by a wetting boundary condition [cf. 
zoom in Fig.[^a)]. Remarkably, we find that the contact 
exhibits at the edge a universal self-similar structure, de¬ 
scribed by a single similarity solution. The similarity 
theory offers a complete description of the physics inside 



FIG. 1: Adhesive contact between an elastic layer and a rigid 
indenter. (a) Stiff contact, the elasto-capillary length 7/p is 
much smaller than the contact size 1. (b) Soft contact, 7/p 
is comparable to 1. Profiles h{x) are solutions of our local 
contact theory (vertical scale stretched). The upper case is in 
the classical JKR regime, for which we identify a narrow, uni¬ 
versal zone that is governed by a wetting angle 8, regularizing 
the crack singularity (inset). 

the contact, including the crack, reconciling the global 
adhesion laws and the local contact mechanics. 

Variational analysis: the wetting angle .— The com¬ 
mon treatment of adhesive contacts hinges on the energy 
released upon opening a crack [in m [Mils]- Here we 
take a different approach by posing the macroscopic free 
energy [35] of the problem sketched in Fig. and derive 
a boundary condition by direct variational analysis. This 
yields a wetting angle 9, instead of the crack singularity. 

As soft materials are nearly incompressible, we need 
to consider only the normal displacement h{x) of the free 
surface of the elastic layer [niiisi- The free energy of the 
problem then consists of an elastic functional, [^] > and 
capillary contributions due to the surface energies 7 , 7 sv, 
and 'ysE [cf. Fig.[^a)]. The associated work of adhesion 
reads W = ^ + 'jsv — Use, representing the reduction 
in the surface energy when contact is created. The free 
energy needs to be minimized under the constraint that 

















2 


the elastic surface must comply with the shape of the 
indenter f{x). This gives the constraint h{x) = f{x) + A 
over the range —i<x<i, where —A is the distance by 
which the solid is indented. Taking also into account the 


work done by the external (two-dimensional) load / 2 D, 
the functional to be minimized is (using the symmetry 
X —>■ —x), 
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The integrals in the upper line give the surface ar¬ 
eas multiplied by their interfacial energies. Compliance 
in the contact zone is imposed by the continuous La¬ 
grange multiplier p{x). Without imposing the slope out¬ 
side the contact, we require h{x) to be continuous at 
a; = £, which is why we introduce the Lagrange mul¬ 
tiplier A to ensure compliance just outside the contact 
[h+{£) = lim 2 ,_,,^+ h(x)]. The resulting variational prob¬ 
lem resembles that of a liquid drop on an elastic layer m- 
The degrees of freedom of the problem thus are 
h{x),£,A, which should minimize the functional Q. As 
shown in the Supplementary Information [28) . the varia¬ 
tions define the elastic problem as [eliminating p{x), A]: 


Cylindrical Indenter .— Within linear elasticity, the 
stress and displacement of the free surface of an incom¬ 
pressible thick layer are related by the integral pS] . 
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valid for all x. We now focus on a cylindrical indenter of 
radius i?, approximated by f{x) = Strict validity of 
linear elasticity requires small strains, i.e. ^ 1 , and 
thus requires that 9 be close to tt. This simplifies the 
boundary condition to h'_ — h'_^ = showing 

that W enters in the form a slope discontinuity. It is 
advantageous to introduce dimensionless variables: 


h{x) 

= /(a;) + A, 

kl < £, 

(2) 

a{x) 

-fh" 

|a:| > £, 
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(1-f /i'2)3/2’ 
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= Isv - ISE, 

X = ±£, 

(4) 


relating to the total force as 


X = x/£, H = hRl£^, F 2 d = ^^, (7) 

TriJ.£^ 

where the external load was scaled by the result for non¬ 
adhesive (Hertz) contacts. It turns out convenient to 
define the following dimensionless numbers: 


ho = 1/2 ~ / dxa{x). (5) 

Here a = SiFeiIdh is the elastic normal stress, which sup¬ 
ports the capillary pressure outside the contact ([^. 

The crux of the analysis is the appearance of Young’s 
wetting angle 9, defined in Fig.[2 it serves as the bound¬ 
ary condition emerging from the variation of £. Like in 
elastic wetting [2TJ Ea [27], surface tension thus domi¬ 
nates over elasticity at the contact line and replaces the 
crack by a “wedge” geometry. Such a wetting-like bound¬ 
ary condition is also hypothesized by Jensen et al. [25] . 
who experimentally found that a gel layer approaches the 
indenter at a well-defined angle that is independent of the 
global geometry. Importantly, from a theoretical point 
of view, the adhesive properties are now contained in the 
theory, since Q can be expressed as IT = 7(1 -I- cos 9). 
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( 8 ) 


S represents the influence of surface tension, comparing 
the elasto-capillary length 7 //J, to the width £ of the con¬ 
tact. A represents the work of adhesion, dictating the 
wetting angle. Following m, we note that outside the 
contact cr(a;) = 'yh"{x), and one can reduce Q to ^5) : 
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valid for A > 1, outside the contact. 

The integral equation ([^ completely defines the ad¬ 
hesive contact problem. Surface tension S is represented 
explicitly, while the adhesion parameter A appears as the 
slope discontinuity boundary condition, H'{\) = 1 — A. 
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FIG. 2: Numerical solutions to the local contact mechanics 
equation § without external load (/ 2 D = 0). The contact 
width t (made dimensionless using the scaling of A) is shown 
as a function of the surface tension 7 (made dimensionless 
using the scaling of S). For small ylfil, we recover the JKR 
scaling [solid line, ( 111 ], while for large 'y/fJ.l the width satu¬ 
rates. The inset shows the corresponding profiles. 


The examples in Fig. [l] are actual numerical solutions 
with S ^ 0. Importantly, the boundary condition cannot 
be imposed without surface tension, 5 = 0, for which the 
solution exhibits the crack-singularity around X = 1, 


i5'(X) ~X-it:(X-l)-i/2. (10) 

Here K is essentially the (scaled) stress intensity factor, 
which for the cylinder reads K = {1 — [26] . 

Results for vanishing load. — To illustrate that our 
local theory indeed provides a complete description, we 
present numerical solutions to with boundary condi¬ 
tion H'{1) = 1 — A, for the case of vanishing load. Fig¬ 
ure]^ shows how the width £ varies with surface tension 7 
[both scaled according to A and S, cf. (|^]. For small S 
we perfectly recover the result for “solid adhesion” (solid 
line) [TTl 113] , which in the present notation reads 


£ = 


8R^W 

TTfJ. 


1/3 


( 11 ) 


By contrast, for large surface tension (or equivalently 
vanishing /t), the contact width becomes independent 
of S and the curve saturates to the value expected for 
“liquid wetting”. This is in agreement with the analysis 
based on energy release rate by Liu et al. [m, but here 
derived from an analysis of the entire contact - including 
the crack region and the angle 9. The result in Fig. 
also concurs with experiments m and simulations m- 
The present theory has the additional merit that it re¬ 
veals the stresses near the contact edge. The inset of 
Fig. [3] shows a stress profile on the scale of the contact 
(5 = 1, black line). The stress extends both inside and 
outside the contact and exhibits a weak, logarithmic sin¬ 
gularity associated to the slope discontinuity [37]. This 



xl■C-^ 


FIG. 3: Stress profiles a outside the contact with £20 = 0, 
for S = lO”"^, 10“^, 10“^, 10“^, 1. The stress exhibits a -3/2 
scaling (dashed line), regularized by a logarithmic zone as 
X ^ 1. Inset: Stress profile on the scale of the contact. Grey 
line is JKR profile (5 = 0), black line a typical profile with 
surface tension (5 = 1). 

is very different from the JKR theory (5 = 0), which 
has the square root singularity at the inside and van¬ 
ishing stress outside the contact (inset, grey line). The 
main panel in Fig. shows further details of the stress 
when approaching x ^ I from the outside of the contact. 
The various curves correspond to different 5. One ob¬ 
serves the appearance of a -3/2 power law, which is loga¬ 
rithmically smoothened at small distances where surface 
tension becomes dominant. The various profiles all look 
very similar. The cross-over observed in Fig. can thus 
be explained from the physics inside the contact: the 
transition from “adhesion” to “wetting” appears when 
the influence of surface tension, the smoothing zone, be¬ 
comes comparable to the size of the contact, i.e. p. ^ 1. 

Universality of the edge region: similarity analysis .— 
We now expose the self-similar nature of the edge region. 
Since the extent of surface tension is bounded to a thin 
region that scales with 5, we introduce a boundary layer 
to regularize the singular crack solution [55]. We there¬ 
fore propose the similarity form, valid for 5^1, 


H{X)=ASn{(:) + lx\ C=^^- (12) 

Here ’H(C) is a universal function describing the inner 
region of the crack. As detailed in the Supplementary 
[28] . inserting ( 12 ) in ([^ gives a closed equation for 'H{Q: 


n'fo = - / 
^0 


dQ 
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(13) 


complemented by the boundary condition 'H'(O) = —1. 

The numerical solution to (13) is represented in Fig. 
showing %"{£,) as the red dashed line. Owing to the capil¬ 
lary stress relation, a = 7 / 1 " = ^AT-L"/RS, the similarity 
Ansatz in fact predicts a collapse of all stress profiles: 
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FIG. 4: The similarity function T-L'^C,) describing the crack 
region (red, dashed). Supe rim posed are the stress profiles of 
Fig. in scaled according to (|14[) , confirming the self-similarity. 




(14) 


Figure]^ shows that the scaled stress profiles of Fig. [^per¬ 
fectly collapse on the predicted dashed line, illustrating 
the validity of the similarity theory. 

Adhesion laws: matched asymptotics .— To complete 
the analysis, we demonstrate how the adhesion laws are 
recovered. We observe that %' ~ —l/(7rC^^^) for large 
C, and this behavior can be matched to the square root 
divergence of (10). This is illustrated in Fig. [^ when 
the effect of surface tension is small. S' <C 1, one observes 
an overlap between the inner solution , red dashed) 
and the outer solution (i/g, gray dashed). Equating the 
prefactors of the asymptotes gives, 




, f2DR 


where in the second step we used the dimensionless stress 
intensity factor K = (I — . 

The matching condition © coincides with the adhe¬ 
sion law for the cylindrical indenter, previously derived 
from a balance of the energy release rate and W mils]. 
For the case jiu = 0 it reduces to (11), while without 


adhesion it gives the Hertz law between the size of the 
contact I and the applied load. The similarity theory 
gives a physical description of this result: the adhesion 
imposes a wetting angle at the edge of the contact, which 
by means of a boundary layer is connected to the elastic 
displacements on the much larger scale of the contact. 

Analogously, the same matching gives the JKR-law for 
the spherical indenter. The boundary layer is asymptoti¬ 
cally thin compared to the contact radius, so the physics 
at the edge is quasi one-dimensional (see also [H]). The 
spherical contact is therefore governed by the same uni¬ 
versal form 'H(C) and again requires a matching to the 


FIG. 5: Matching the asymptotes of the similarity solu¬ 
tion %' (red, dashed), to the singular outer solutions H'^/A 
for the cylinder (gray dashed) or sphere (blue dotted). The 
solid black line shows the 2D numerics: it follows Hq on large 
scales, switching to H' on small scales. The matching was 
done for S = 10“^^, while the external load was set to the 
critical value for pull-off (inset, circle). Inset: £ versus fsn 
given by Eq. (161, scaled by the critical values fc,£c. 


singular Hq (Fig. [^ blue dotted). Evaluating K for the 
spherical indenter indeed gives the JKR-law [5]: 


K = 


ihoR /97rWR2y/2 

16/r£3 Iv ) 


(16) 

Eor completeness, we illustrate this famous relation be¬ 
tween £ and the load in the inset of Eig.[^ The outer 
solutions shown in the main panel of Eig. [^were selected 
to correspond to the case where the external force reaches 
its minimum (negative) value, which is the critical point 
where the indenter is pulled off from the elastic layer. 


Outlook. — The present analysis provides a complete 
unification of solid adhesion and wetting: both are gov¬ 
erned by a wetting angle, and this regularizes the crack 
singularity. While the elastic solution given here is re¬ 
stricted to linear elasticity, (VF/y)^/^ <C 1, the boundary 
condition (j^ equally applies nonlinearly - as also sug¬ 
gested by recent experiments [25] . A noteworthy nonlin¬ 
ear effect at the contact line is that the solvent can be 
separated from the gel in the vicinity of the edge , in 
particular close to complete wetting [50] • In a broader 
context, we hypothesize that self-similarity could be a 
generic feature of fracture in the presence of surface ten¬ 
sion ISTl |32| , as is the case for fracture of viscous liquids 
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